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ABSTRACT 

We present a state-of-the-art linear redshift distortion analysis of the recently pub- 
lished IRAS Point Source Catalog Redshift Survey (PSCz). The procedure involves 
linear compression into 4096 Karhunen-Loeve (signal-to-noise) modes culled from a 
potential pool of 3 x lO'^ modes, followed by quadratic compression into three sep- 
arate power spectra, the galaxy-galaxy, galaxy-velocity, and velocity-velocity power 
spectra. Least squares fitting to the decorrelated power spectra yields a linear redshift 
distortion parameter (i « ^^m^/f* = 0.411 
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1 INTRODUCTION 

The IRAS Point Source Catalog Redshift Survey (PSCz), 
made public on 9 January 2000 (Saunders et al. 2000), cov- 
ers more volume of the Universe than any other publically 
available redshift survey. This, together with the consider- 
able care taken by its authors to ensure uniformity of selec- 
tion over a prescribed volume, makes the PSCz the finest 
available galaxy survey for measuring the power spectrum 
and its redshift distortions at large, linear scales. 

The purpose of this Letter is to report a state-of-the- 
art finear redshift distortion analysis of the PSCz. Our ap- 
proach has its roots in the work of Fisher, Scharf & Lahav 
(1994), who were the first to apply a maximum likelihood 
approach to the analysis of linear redshift distortions, and 
in the work of Heavens & Taylor (1995, hereafter HT), who 
may be credited with accomplishing the first linear likeli- 
hood analysis designed to retain as much information as 
possible at the largest, linear scales. Fisher et al. and HT ap- 
plied their approaches to the the IRAS 1.2 Jy survey (Fisher 
et al. 1995), obtaining values of the redshift distortion pa- 
rameter /3 = 0.96jig jg and j3 — 1.0 ±0.3 respectively, where 
f} = //6 is the dimensionless linear growth rate / ~ Q,'^^ 
divided by the linear galaxy-to-mass bias factor b. A com- 
prehensive review of these and other measurements of linear 
redshift distortions, complete up to mid-1997, may be found 
in Hamilton (1998). 

Most recently, Tadros et al. (1999) applied HT's method 
to the PSCz survey, conservatively cut to 0.75 Jy. To deal 
with the more complicated PSCz mask, Tadros et al. used an 
improved treatment of the angular mask, allowing them to 
include monopole and dipole modes. Tadros et al. obtained 
(3 — 0.58 ± 0.26 (marginalized error) if the shape of the 



power spectrum was fixed to be that of a F = 0.2 CDM 
model, or /3 = 0.47 ±0.16 (conditional error, with the power 
spectrum fixed at its best fit) if the power spectrum was 
simultaneously measured in 6 bins. 

The goal of the present paper is to measure the linear 
power spectrum and its redshift distortions at high reso- 
lution with minimal loss of information. The paper brings 
together several ideas from the literature, chief amongst 
which are: (1) use of the logarithmic spherical wave ba- 
sis (Hamilton & Culhane 1996); (2) linear compression into 
Karhunen-Loeve (KL) , or signal-to- noise, eigenmodes (Bond 
1995; Bunn & Sugiyama 1995; Vogeley & Szalay 1996; Mat- 
subara, Szalay & Landy 2000); (3) quadratic compression 
into high resolution powers (Tegmark 1997, 1998, hereafter 
T97, T98; Tegmark, Taylor & Heavens 1997, hereafter TTH; 
Tegmark et al. 1998; Padmanabhan, Tegmark & Hamilton 
2000); (4) inclusion of stochastic bias in the analysis (Pen 
1998; Tegmark & Peebles 1998; Dekel & Lahav 1999); (5) 
decorrelation (Hamilton 1997; Tegmark & Hamilton 1998; 
Hamilton & Tegmark 2000); (6) the elimination of bias re- 
sulting from the 'pair-integral constraint', by the trick of 
making all KL modes orthogonal to the mean mode (Fisher 
et al. 1993; Tegmark et al. 1998). 



2 THE 2-POINT CORRELATION FUNCTION 
IN REDSHIFT SPACE 

Figure 1 shows the redshift space 2-point correlation func- 
tion of PSCz. The expected large scale squashing effect 
caused by infall toward overdense regions is plainly visible. 
Nonlinear fingers-of-god show up as ~ 4/i~^Mpc wide by 
20/i~^Mpc long enhancements along the line-of-sight axis. 
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Figure 1. Contour plots of the redshift space two-point corre- 
lation function in the PSCz survey with the high galactic lati- 
tude angular mask, (left) beyond 25/i~^Mpc, and (right) beyond 
80/i~^Mpc. Thin, medium, and thick contours represent nega- 
tive, positive, and zero values respectively. The correlation func- 
tion has been smoothed over pair separation r = (r^ + rp'^/^ 
with a tophat window of width 0.2 dex, and over angles 9 = 
tan~^ (rx/^ii ) to the line of sight with a Gaussian window with 
a Icr width of 10°. 

There is little sign of the large scale 'arm-of-god' af- 
fliction beyond 80/i~^Mpc found by Hamilton (1995) in the 
QDOT survey, the l-in-6 precursor to PSCz (Lawrence et 
al. 1999). The contour plots in Figure 1 are constructed in 
a manner essentially identical to those in Hamilton (1995) 
and Hamilton (1998). The only difference is that, thanks to 
the lower noise level of the PSCz, the contour plots in Fig- 
ure 1 extend to larger scales and to lower contour levels than 
those of Hamihon (1995, 1998). 



3 ANALYSIS 

The analysis will be described in full in a subsequent paper. 
Here we summarize the main points. 

As emphasized by Tadros et al. (1999), in studies of 
weak clustering at large scales, uniformity of selection is 
paramount, since non-uniformity will masquerade as spu- 
rious excess power. Following the advice of Saunders et al. 
(2000) and Tadros et al., we adopt the high- latitude angular 
mask (hibpsczmask.dat, part of the PSCz package). 

It is important to measure the radial selection function 
of the survey as accurately as possible, since errors in the se- 
lection function translate into spurious large scale power. We 
use the maximum likelihood method of Sandage, Tammann 
& Yahil (1979), and fit the selection function to a smooth 
analytic function with enough free parameters to yield an ex- 
cellent fit. The assumed geometry is that of a flat flm = 0.3, 
SIa ~ 0.7 Universe. The measurement yields evidence for 
what appears to be strong evolution, which we model as lu- 
minosity evolution with luminosity oc (1 + z)^'"^. Similarly 
large evolution was previously reported in the QDOT sur- 
vey by Saunders et al. (1990). Figure 2 shows the observed 
number density of galaxies, divided by the measured selec- 
tion function, both with and without evolution. 

Figure 2 also shows what appears to be incompleteness 
beyond ~ 420/i~^Mpc. This may be presumed to be the 




Figure 2. Ratio of the observed galaxy number density to the 
maximum likelihood selection function at radial depth r in the 
PSCz survey, averaged in depth bins 0.025 dex wide. The solid 
line assumes that galaxies evolve with luminosity oc (1 + z)^'^, 
while the dashed line assumes no luminosity evolution. For the 
analysis of this paper, the survey is cut at a comoving depth of 
;^q2.625 /j-lMpc 420/i~lMpc, indicated by the vertical dashed 
line. 

incompleteness at high redshift described in §4.2 of Saun- 
ders et al. (2000), associated with the policy not to pursue 
redshifts of galaxies optically fainter than bj = 19.5™. Since 
this incompleteness is greater in regions of higher optical ex- 
tinction, and is systematic rather than random over the sky 
(Fig 4. of Saunders et al. 2000), we choose to cut the survey 
at lO^ ''^^ tr^Mpc « 420/i~^Mpc. We also choose to cut the 
survey at a near distance of 1 /i~^Mpc, in order to eliminate 
the local nonlinear region, the Local Group of galaxies. This 
leaves 12504 galaxies in the sample. 

As a working basis from which Karhunen-Loeve modes 
(see below) are constructed, we use logarithmic spherical 
waves (Hamilton & Culhane 1996), which are orthonormal 
eigenfunctions Z^<„(r) = (27r)-i''2e-(3/2+i-)i°'-y^„(^) of 
the complete set of commuting Hermitian operators 

The advantages of this basis are: (1) the linear redshift dis- 
tortion operator is diagonal in this basis (Hamilton & Cul- 
hane 1996); (2) radial modes can be transformed rapidly 
between real, uj, and Fourier space using Fast Fourier Trans- 
forms; (3) the radial modes discretize naturally on to a grid 
of radial depths r that is uniformly spaced in the log, which 
makes it well suited to surveys like PSCz, where the flux 
limit causes the survey to be finely sampled nearby, and 
coarsely sampled far away. 

The principal numerical limitation to the linear likeli- 
hood method is the number TV of modes that can be included 
in the Gaussian likelihood function. Solving the likelihood 
formulae involves manipulating N x N matrices, which is a 
~ A'^^ process. Thus the modes should be crafted so as to 
include as much information as possible about the quantities 
of interest, the parameters to be measured (HT; TTH). 

Here we use Karhunen-Loeve (KL), or signal-to-noise, 
modes, as first suggested in the context of galaxy cluster- 
ing by Vogeley and Szalay (1996). The covariance matrix 
of galaxy densities is a sum C = S + N of cosmic signal S 
and Poisson sampling noise N. KL modes are constructed 
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by diagonalizing the signal-to-noise matrix N-i/2g]sj-i/2 
and selecting the eigenmodes with the largest eigenvalues. 

The main difficulty with the KL procedure is that in 
order to construct N modes it is necessary to diagonalize a 
matrix of dimension ^ N . But this would seem to defeat 
the original goal, to avoid manipulating huge matrices. We 
get around this difficulty by constructing KL modes in a two 
stage process, constructing first angular KL modes, then rar 
dial KL modes within each angular mode. The resulting KL 
modes arc not perfect, since they ignore covarianccs between 
radial modes belonging to different angular modes. But the 
KL modes need not be perfect; it is enough that they should 
contain virtually all the information of interest. 

In practice, we use 1600 angular modes, and 192 radial 
modes, so there is a potential pool of 1600 x 192 ~ 3 x 10^ 
modes. This is 75 « 4? times the number 4096 of KL modes 
retained in the analysis. In eflfect, the KL modes are spatially 
over-resolved by a factor of about 4 in each dimension. 

In order to isolate the effects of the selection function, 
and of the motion of the Local Group (LG) through the Cos- 
mic Microwave Background (CMB), we choose the first four 
angular 'KL' modes to be the cut (i.e. masked) monopolc, 
and three cut dipole modes. The cut dipolc modes contain 
admixtures of the cut monopole, to make them orthogonal 
to the latter, and all other angular KL modes are orthogonal 
to the cut monopole and dipole. Within the cut monopole 
mode, the first two radial modes are the mean (selection 
function) mode and the Local Group mode (equation 4.42 
of Hamilton 1998), the latter containing a small admixture 
of the mean to make it orthogonal to the cut (i.e. radially 
masked) mean. Within the cut dipole modes, the first radial 
mode is the Local Group mode. 

The mean mode is used in determining the maximum 
likelihood normalization of the selection function, but is then 
discarded from the analysis, since it is impossible to measure 
the fluctuation of the mean mode. 

Since all KL modes are orthogonal to the mean, there 
is no need to subtract the mean. Making the KL modes or- 
thogonal to the mean immunizes their amplitudes against 
uncertainty in the mean density, and ensures that the mea- 
surement of power is unbiased by the so-called 'pair-integral 
constraint'. This trick was first used by Fisher et al. (1993). 

The fact that the selection function is measured in red- 
shift space, not real space, leads to a correction to the lin- 
ear redshift distortion operator for the monopole modes, de- 
scribed in §4.4 of Hamilton (1998). The correction leads to a 
noticeable reduction in measured power at the largest scales. 

The motion of the LG through the CMB is known 
(Courteau & van den Bergh 1999; Lineweaver et al. 1996). 
We correct the amplitudes of the four LG modes (one cut 
monopole and three cut dipole) for this motion, and include 
them in the analysis, as did Tadros et al. (1999). 

The amplitudes of the 4096 KL modes constitute the 
'data' which now feed into the Gaussian likelihood func- 
tion. Rather than attempting to solve the likelihood function 
directly for values of cosmological parameters, we instead 
quadratically compress (T97; TTH) the data into an inter- 
mediate number (here 147, plus one for the normalization of 
the selection function) of parameters that contain essentially 
all the information of interest. The advantage of quadratic 
compression is that it allows large numbers of parameters 
to be measured rapidly, without a time-consuming nonlin- 



ear search for an extremum in a high-dimensional space. 
Key to the quadratic method is that the intermediate pa- 
rameters should be linearly related to the mean and covari- 
ance. In a linear redshift distortion analysis, the intermedi- 
ate parameters that emerge naturally are the galaxy-galaxy, 
gala^xy-velocity, and velocity-velocity power spectra (Kolatt 
& Dekel 1997; Tegmark 1998; Pen 1998; Tegmark & Peebles 
1998; Dekel & Lahav 1999): 

galaxy-galaxy power : -Pgg(fc) = b{k)^P{k) 
galaxy-velocity power : -Pgv(fc) = r{k)b{k)fP{k) (2) 
velocity-velocity power : Pvv{k) = f^P{k) 

where b{k) is the (possibly scale-dependent) galaxy-to-niass 
bias factor, r{k) £ [—1,1] is a (possibly scale-dependent) 
galaxy- velocity correlation coefficient, / « is the di- 
mensionless linear growth rate, and P{k) is the mass power 
spectrum. More correctly, the 'velocity' here refers to minus 
the velocity divergence, which in linear theory is related to 
the mass (not galaxy) overdensity S hy fS + V -v = 0, where 
V denotes the comoving gradient in velocity units. 

Although the three power spectra are interesting in 
their own right, here we report results in terms of conven- 
tional parameters: the linear redshift distortion parameter 
(3 = f/b, the galaxy-velocity correlation coefficient r, and 
the galaxy-galaxy power spectrum i-'gg(fe). 

The estimate of power that emerges most directly from 
quadratic compression is the power spectrum convolved with 
the Fisher matrix (T97; TTH). Physically, this represents 
the power spectrum convolved with the Fourier transform 
of the optimally weighted survey window. We call this es- 
timate, appropriately renormalized, the 'correlated' power, 
since the errors are correlated. Deconvolving the correlated 
power yields the 'raw' power spectrum. The raw power tends 
to be noisy and anti-correlated, and is not very useful. Half 
way between the correlated and raw power spectra is the 
'decorrelated' power spectrum, which is partially decon- 
volved in such a way that estimates of power at different 
wavenumbers are uncorrelated with each other (Hamilton 
1997; Hamilton & Tegmark 2000). 

The quadratic method requires that prior powers 
Pgg(fe), Pgv(fc), and Pvv(fc) be specified. The maximum like- 
lihood (ML) solution is that power spectrum for which the 
estimated power equals the prior power. The ML solution 
can be obtained by folding the estimate back into the prior 
and iterating to convergence (Bond, Jaffe & Knox 1998). 

However, in this paper the power is permitted to be an 
arbitrary function of wavenumber, with no penalty against 
violently varying power spectra. To express our Bayesian 
prejudice in favour of a smoothly varying power spectrum, 
we choose to fold back into the prior not the raw power 
spectrum, but rather the correlated power spectrum, which 
is smooth(ish). Moreover we iterate only once, since there is 
no point in attempting to overfit the bumps and wiggles in 
the power. As the initial prior power, we choose a COBE- 
normalized flat ACDM model with Sim ~ 0-3, SIa = 0.7 from 
Eisenstein & Hu (1998). 

Further, we set the redshift distortion parameter in the 
prior power to a constant, /? = 0.5, reflecting our Bayesian 
prejudice that this parameter should be constant in the lin- 
ear regime. The data do not constrain the galaxy-velocity 
correlation coefficient r tightly, but are consistent with r « 1 
at linear scales, so we fix r = 1 in the prior. 
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Figure 4. Linear redshift distortion parameter fS = f/b as a function of wavenumber k in the PSCz 0.6 Jy survey. (Left) Dashed 
line is estimated rf3 from the ratio of galaxy- velocity to galaxy-galaxy power, Pgv/Pgg. Solid line is estimated /3 from the square root 
of the ratio of velocity-velocity to galaxy-galaxy power, sign(Pvv/^'gg)|Pvv/^gg|^''^- (Middle) Best fit value of /3 from the combined 
ratios of galaxy-velocity and velocity-velocity to galaxy-galaxy powers, assuming unit galaxy-velocity correlation coefficient, r = 1. For 
comparison, the graph also shows the value of /3 inferred from the ratio of quadrupole-to-monopole power, computed in the manner of 
Hamilton (1995, 1998). (Right) Cumulative best fit value of fS including all data at wavenumbers less than or equal to the wavenumber 
specified on the x-axis, again assuming r = 1. The shaded region is the l-cr uncertainty in /3, marginalized over the overall amplitude 
of the power spectrum. 
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Figure 3. Galaxy-galaxy power spectrum measured from the 
PSCz 0.6 Jy survey with the high latitude angular mask. This 
is the real space power spectrum, the redshift distortions hav- 
ing been isolated into galaxy- velocity and velocity- velocity power 
spectra. The solid line is the correlated power spectrum, and the 
shaded region its l-cr error, while points with error bars constitute 
the decorrelated power spectrum (Hamilton &; Tegmark 2000). 
Each point of the decorrelated power spectrum is uncorrelatcd 
with all other points. The dashed line is a linear model power 
spectrum from Eisenstein & Hu (1998), a COBE-normalized, un- 
tilted, flat ACDM model with Cm = 0.3, Qa = 0.7. 



4 RESULTS 

Figure 3 shows the galaxy-galaxy power spectrum measured 
from the PSCz 0.6 Jy survey. This is the real space power 
spectrum, the redshift distortions having been isolated into 
the galaxy-velocity and velocity-velocity power spectra. At 
nonlinear scales, k ^ 0.3 feMpc"^, the interpretation as a 
real space power spectrum becomes suspect, since the extent 
to which nonlinear redshift distortions can be approximated 
as (fc-dependent) linear combinations of the linear galaxy- 
velocity and velocity-velocity powers is unknown. 

Figure 3 also shows a linear model power spectrum 



from Eisenstein & Hu (1998), an observationally concor- 
dant COBE-normalized, untilted, flat ACDM model with 
f2m ~ 0.3, f^A ~ 0.7, baryonic content Qf,h^ = 0.02, 
and Hubble constant h — 0.65. Interestingly, the ACDM 
power spectrum fits the PSCz data well at linear scales, 
k 0.3/iMpc~^, with no bias, b — 1. The enhancement of 
the observed power over the model power at smaller scales, 
k J; 0.5 hMpc~^, can be attributed to nonlinearity (Peacock 
& Dodds 1996). 

The real space power spectrum measured here appears 
consistent with the redshift space power spectrum measured 
from PSCz by Sutherland et al. (1999). 

Figure 4 shows the redshift distortion parameter (3 = 
f/b measured as a function of wavenumber k. The left panel 
shows the values of (3 measured separately from the ratio 
Pgv(fe)/Pgg(fe) of the galaxy- velocity to galaxy-galaxy power 
on the assumption of unit galaxy-velocity correlation coeffi- 
cient, r = 1, and from the ratio sign(Pvv/fgg) \Pvv / Pgg\^'''^ 
of the square root of velocity-velocity to galaxy-galaxy 
power, equations (2). These are ratios of correlated power 
spectra, since the decorrelated galaxy-velocity and velocity- 
velocity powers are confusingly noisy. 

The central panel of Figure 4 shows the best fit value 
of /3 that comes from combining the two ratios of galaxy- 
velocity and velocity-velocity to galaxy-galaxy power at each 
wavenumber k. Also shown for comparison is the value of 
/3 inferred from the ratio of smoothed quadrupole power 
to smoothed monopole power, measured in the fashion of 
Hamilton (1995, 1998). The two measures of /3 agree at scales 
k ^ 0.04/iMpc"\ The difference at fc ^ 0.04 /iMpc"^ may 
plausibly be attributed to the inability of the quadupole-to- 
monopole ratio to probe the largest scales. 

The right panel of Figure 4 shows the best fit value of 
j3 that emerges from a least squares fit to the decorrelated 
powers including data cumulatively up to the wavenumber 
specified on the a:-axis, again on the assumption of unit 
galaxy- velocity correlation coefficient, r = 1. The shaded 
region bounds the l-cr uncertainty in /3, where exceeds 
its minimum value by 1. 
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In carrying out the fit we fix tlic sliapc of tlic galaxy- 
galaxy power spectrum to be that of the measured correlated 
power spectrum shown in Figure 3, but we allow the overall 
amplitude of the power spectrum to vary. The l-cr uncer- 
tainty in Figure 4 is the uncertainty marginalized over the 
amplitude of the power spectrum. We also tried allowing the 
shape of the power spectrum to vary, but the shape is al- 
ready measured about as well as can be, and we could not 
thereby reduce in a statistically significant way (reduc- 
tion by at least 1 for each additional degree of freedom) . 

The fit is consistent with constant (3 tor k ^ 
0.3 ?iMpc~^, and as an overall best fit we quote the value 
l3 = 0.4lj:g:i^ at fc = 0.274 /iMpc-^ ACDM with Q^n = 0.3 
and = 0.7 predicts / = 0.513, consistent with our best 
fit result if IRAS galaxies are unbiased or mildly biased. 

We also tried allowing the galaxy-velocity correlation 
coefficient r to vary. At linear scales, k ^ 0.15 /i Mpc"'^, the 
best fit value is of order unity, r « 1, as suggested by the 
approximate agreement in the two curves in the left panel 
of Figure 4, but the uncertainty is largo. At smaller scales, 
k ^ 0.15/iMpc~^, the measured coefEcient decreases be- 
low 1. This is almost certainly caused by nonlinearity, not 
by stocliasticity in the linear galaxy-velocity correlation. In 
fact fingers-of-god, which appear extended along the line-of- 
sight in redshift space rather than compressed, are expected 
to cause the galaxy-velocity power to go negative at nonlin- 
ear scales, yielding a negative galaxy-velocity correlation co- 
efficient r. The left panel of Figure 4 shows that the galaxy- 
velocity power indeed goes negative at k ^ 0.2/iMpc~'^. 



5 CONCLUSIONS 

This Letter reports a measurement of the linear power spec- 
trum and its redshift distortions in the IRAS PSCz 0.6 Jy 

survey (Saunders et al. 2000) using state-of-the-art meth- 
ods, including logarithmic wave expansion, Karhunen-Loeve 
compression, quadratic compression, decorrelation, and iso- 
lation of the mean mode in order to eliminate large scale 
bicis from the 'pair-integral-constraint'. 

The measurement yields three separate linear power 
spectra, the galaxy-galaxy, galaxy-velocity, and velocity- 
velocity power spectra. 

The best fit overall value of the linear redshift distortion 
parameter, /3 = f/b = 0.41^q;j2) is consistent with, albeit 
slightly lower than, that reported by Tadros et al. (1999) for 
a more conservative 0.75 Jy cut of the survey. 

The galaxy-velocity correlation coefficient r is not 
tightly constrained by the data, but is consistent with being 
about unity at linear scales, r « 1. 

The inferred real space power spectrum and its red- 
shift distortions are consistent with a COBE-normalized, 
untilted, flat ACDM model with = 0.3 and Qa = 0.7, if 
IRAS galaxies are an unbiased tracer of mass at large scales 
in the Universe. 
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